DAMPING ROTOR NUTATION OSCILLATIONS 
IN A GYROSCOPE WITH MAGNETIC SUSPENSION 1 ' '' u 


Valentine N. Komarov 

Research Institute for Applied Mathematics and Cybernetics 
of Nizhny Novgorod State University, Nizhny Novgorod, Russia 


SUMMARY 


A possibility of an effective damping of rotor nutations by modulating the field of the 
moment transducers in synchronism with the nutation frequency is considered. The algorithms 
for forming the control moments are proposed and their application is discussed. 


INTRODUCTION 


Practical use of noncontact suspensions in gyroscopy is possible if the following 
requirements, at least, are fulfilled: the rotor must be stably suspended and stably rotate with 
given parameters; information about the angular motion of the body must be available. 

The second problem of making the rotor rotate is, in turn, divided into three other problems: 
providing a necessary angular momentum for the rotor, which, most frequently, is equivalent to 
providing a certain rotation speed Q for the rotor, providing a small, zero in the limit, nutation 
angle d , and maintaining the achieved motion parameters within the given limits in the further 
operation of the unit. A way for solving these problems with respect to a gyro with a 
magnetoresonance suspension can be found in (ref. 1). 

The small dissipative moments in a noncontact gyro leads to the fact that after acceleration 
the rotor rotates with a nutation angle defined mainly by the initial conditions. The damping of 
the nutation oscillations over a finite time can be achieved if, for the formation of the control 
moments (not only for nutation damping), we use modulation or self-modulation methods, when 
parameters of the suspension field or of the additional coil field (magnitude, direction, rotation 
speed, etc.) vary depending on an angular position of the aspheric rotor (ref. 2, 3, 4), due to 
which a nonconservative component appears in the moment effect on the rotor, which also 
provides an active damping of nutation. Let us consider one of the possible methods for an 
effective damping of the rotor nutation relative to a gyro with magnetic suspension. 


‘This research has been done under the support of the Russian Foundation for Fundamental 
Research (Project 95-01-00002) 
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THE EQUATIONS OF MOTION 


Assuming the gyro rotor to be stably suspended in the suspension field, we represent the 
equations of its angular motions in the form of six first order equations with respect to the 
following phase variables (ref. 5): the value of angular momentum I\ , two angles p and a . 
wdiich determine its orientation (trihedron OYj) relative to the suspension (trihedron OZ\ t), and 
The Euler angles p , i p and d , which prescribe the position of the rotor (trihedron (O.V,)) 
relative to the angular momentum K ; p is the angle of the proper rotation of the rotor about 
the axis of dynamic symmetry {OX3) ; tp is the angle of precession of this axis about the 
angular momentum K j i? is the nutation angle or the angle between the angular momentum and 
the axis OX 3. 

Using the main equation of gyroscopy 

dK 

— |fixK = M ( 01 ) 

written in the coordinate system OYj and projecting the vector K on the rotor axes OX, 

(K • X.) = K • 6,3 = I, ft, , (02) 

where /, is the rotor’s moment of inertia along OXi, and ft; is the projection of the rotor’s 
absolute angular velocity on the same axis, for an axisymmetric rotor (/, = J 2 = A ^ / 3 = C) 
and the arbitrary moments we obtain 

K ( J + ^’ 2 ) = ’ K ( sin/> ^ - A.i j = M 2 , ~ = M 3 , (03) 

,dd 

A — = - (A/, cos V’ + M 2 sin V’) = -M con , 


dip 

dt 



1 \ n M 2 cos V’ — Mi sin xb 

— I cos v = — 

A) h sin 1? 


dij> A A/ 2 cos ip — M t sin V’ 
~dt ~A ~ K 


cotan d — ft 3 , 


where Mj — • y , ft 3 — ft.,3 + ( M 2 / A ) cotan/7 , and ft ,j are the piojcct ions onto the 


axes 0 } j of the absolute angular velocity ft, of the coordinate system OZ This velocity can be 
determined by the earth’s rotation, by the rotation of an object on which the gyro is installed, 
by a forced rotation of the device casing performed to make its interaction with the rotor 
symmetric, etc. 


The presence of different time scale motions in real gyros, such as a "slow” precession motion 
of the angular momentum and a ’’fast” rotor rotation allows 11s to use, for studying the gyro 
dynamics, the averaging method in which even the first approximation reveals the basic laws of 
gyro behavior (ref. 5). 
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Omitting the standard procedure for normalizing the equations and finding in the latter a 
small parameter £,. which characterizes the ratio of the work of external forces per one rotation 
of the rotor to the kinetic energy accumulated by the gyro rotor, we write their standard form, 
which allows the asymptotic methods to be applied: 

p* = e , <r* sin/> = e , k* = Eiiij , (04 ) 

kfl* = —e (rnj cos V’ + m 2 sin iji) , 


<p* -f 7 k cos i? = e 


' m 2 cos 0 — mi sin </’ 


P - k = -e ^ 


fc sin 

\ 

m 2 cos 0 — 77i i sin if? 


K 


cotan d 


— U)3^ . 


It is evident from the equations that free motion of the rotor (generating solution) represents 
a regular precession or the Euler-Poinsot motion with the constant angular velocities of the 
proper rotation ip* and precession for a constant nutation angle r), such that 


V’o — > ^0 — cos i?o i 



C-A\ 

A ) ’ 


(05) 


where k 0 is the constant value of the angular momentum k which is immobile in the inertial 
space. For the periodic dependence of the right-hand sides of the equations of motion on the fast 
variables xji and p , to obtain the first approximation over the small parameter e , the right-hand 
sides of the equations for the slow variables are averaged over the variables V’ and p , whose 
dependences on r are assumed to be the same as those in the unperturbed motion. 


Let us now specify the requirements to be met by the moment controlling the rotor motion 
to perform the damping of the rotor’s nutation oscillations over a finite period of time. 


After the averaging of the equations of motion over the fast, variables i/' and p , the nutation 
angle varies according to 

kd* — — £m C on (/>, <r, , k) . (06) 


Since the initial values of the nutation angle d are closer to 0 rather than tt/ 2 , and this 
angle must still decrease with time, we introduce variable x = sin t? and linearize Eq. (06) with 
respect to this variable 


kx * = — e 


m°(p,<7, i?, k) + <r,i?,fc)j 


where = m co „{p, a, 0, k), m l = dw con /dx\ r=0 . 


Therefore, a monotonic decrease of the nutation angle to zero over a finite period of time is 
possible only if the condition m° > It?) 1 | is fulfilled. In this case, there is no equilibrium state for 
the angle 6 or r , and their variation rates are always negative, which allows the required effect 
to be achieved. 
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THE FORMATION OF THE CONTROL MOMENT 


Let us discuss a possibility for damping the rotor nutations through its interaction with the 
field of a transducer for the angular position of the rotor. In a gyro with a magnetoresonance 
suspension, this transducer is embodied as a system of two four leg magnetic circuits located on 
the same axis but on the different sides of the rotor with a coil (ref. 1 ) at each leg. An axial 
opening in the rotor leads to the fact that from the moment standpoint, the rotor will behave as 
aiaspheric body, whose surface is described by a set of even harmonics (even Legendre 
polynomials), the first of whicli is the second harmonic responsible for the ellipsoidality or the 
uniaxial anisotropy. The moment defined by this harmonic will be considered below. 

This moment is described by 

M = - 2a (s ■ h) [s x h] , (07) 

where a > 0 is the interaction amplitude determined, in particular, by the field amplitude: s and 
h are the unit vectors of the symmetry axes of the ellipsoid and the field. (Since we allow for the 
anisotropic properties of the rotor determined by its opening, from the magnetic viewpoint it 
should be considered as an oblate ellipsoid tending to position itself with its axis perpendicular 
to the field’s axis, which is reflected by the minus sign in Eq. (07).) 

It can be assumed that the unit vector s coincides with the unit vector of the dynamic 
symmetry axis OX 3 of the rotor, and, thus, 

s = X 3 = yi sin d cos 4’ + y2 sin d sin 4 > + y 3 cos d . 


Using the spherical angles a n and , we prescribe the position of the unit vector h . which 
determines the symmetry axis of the magnetic field generated by the n-th pair of the coils 
belonging to different magnetic circuits but located on the same axis that passes through the 
renter of the rotor 

h = z t sin n„ cos /?„ + z 2 sin a,, sin U + z n cos e» n , 
where o„ = a is the same for all n and f) n = /? 0 + trn/ 2 . 

Let us substitute these unit vectors into the expression for the moment (07) and find its 
component m con = cost/’ + rn 2 sin ’/’ responsible for the variation of the nutation angle d 

M ccn = 2a n (hi sin 4' — h 2 cosy’) A 3 cos d— (08) 

— | |‘ 2 A| h 2 cos2y!’ + {h\ - A,) sin2y’j sin 2d , 

where 

hj = (h ■ }0 , A 1 = cos p sin a cos (rr — j3 n ) — sin pros a . 

A 2 = - sin a sin (a - /?„) , A 3 = sin psin o cos (rr — 3 n ) + cos pros a . 

From Eq. (08) it follows that \f. on is a periodic function of the precession angle »> without 
the constant component A/ c ° on (p,<r, 0 , /l) , which can appear only in the case of a pulsed 
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switching-on of the field over the time periods, during which the moment does not change its . 
sign. As it is evident from Eq. (08), the moment has a harmonic component with frequency V 
proportional to cos 2 t? and a component with double frequency 24>* proportional to sin2i?. For 
i? = 0 this results in 

M con = 2a„ (hi sin xp - h 2 cos ip) h 3 . 


Apparently, for the effective damping of nutation, i.e., dropping the angle d down to zero 
over a finite period of time, the coil field with frequency V’’ must be changed. (The speed of the 
rotor rotation about the symmetry axis OX 3 is equal to Jlo = K/C , the nutation frequency 
V’* = K/A , and, consequently, ip* = C/A ). It should be noted that nutation oscillations of a 
cryogenic gyro can be damped in a similar manner (ref. 3). 

Let us write hi sin ip — h 2 cos tp in the form: 

hi s'mip — h 2 cosxl’ = — sin cos a sin 4’ + 

+ i [(1 + cos p) sin (ip + <T - Pn) - (1 - cos p) sin (i' - cr + /?„)]sin« . 


From the above expression it follows that the optimal switching-011 of the field should follow 
the algorithm 

a, 2Jbr < i’ + a-Pn < (M + I)*- 


a(4’) = 


(09) 


0, (2ik + 1 )tt < V’ + ° ~ Pn < 2(A- + l)ir , 
i.e., the field is switched on over half the period of the nutation oscillations with frequency xp . 


For this operation mode, two pairs of the control coils are switched on at any instant, which 
greatly improves the efficiency but leading, however, to a higher heat release that can be reduced 
due to a decrease in either the field amplitude, or the switching-on duration. 


a(i/’) = 


( 10 ) 


a, (2k +|)7 T-|<V’ + °'-^n< (2fc + |) n + 2 
f 0, outside the intervals 

Here r < w is the switching-on duration. Apparently, the algorithm (09) is a particular case of 
( 10) for r = 7T. 

If all the coil pairs and their switching-on algorithms are identical, they create a total control 

moment T 

(Mcon) = 2- ( 1 + cos p) (2 cos p - 1 ) sin 2a sin - cos d , (11) 

7T * 

which leads to equation 

K 'd* = — 2 — (1 + cos p ) (2 cos p — 1) sin 2ot sin — cos 2 1? . (12) 
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Apart from the variable d , this equation also incorporates other quantities A (f) and p{t) 
that do not remain constant. If we write the rotor- affecting moments via the force function II' 
(ref. 5), then, instead of the corresponding equations in (03), we can obtain 

..dp 1 (dW 5tV\ 

A dt ~ sin p V 50 C ° Sf> ~ ’ 

d K dW 
dt dtp ' 

dd i fdw 5ir\ 
dt K sin d y dip dip J 


Since the rotor opening is symmetric about its dynamic axis, then, consequently, dW/dip = 0 
. The summation of the moments from all the pairs of the symmetric coils also results in 
dW/dcr = 0 . These conditions yield the equations 

..dp 1 dW dl< dW dd 1 d\V 

dt sin p dip COS P ’ dt 50 ’ ’ dt I\ sin d 50 

from which we obtain two integrals 


Kz3 = K cos p = Const , A r ) = K cost? = Const , 


(13) 


the physical meaning of which is evident, i.e., on the strength of the control system and rotor 
symmetries, the projections of the angular momentum onto the symmetry axes of the suspension 
( A 0 t ) and the rotor ( A 03 ) remain constant. Representing K and cos p from Eq.(13) and 
substituting them into Eq. (12), we can find the law of the nutation angle d variation with time. 
To simplify the solution, let us use the smallness of the misalignment angle p , which, along with 
the use of the second integral, yields 


d* = — v cos 3 d , 


4 a . . r \ 

v = — — — sin 2a sin — 

ttA r3 2 ) 


( 11 ) 


from which it is evident that the time of the nutation angle decrease from the initial value d 0 to 
0 is given by 


T(d 0 ) 


1^ f sin dp Jl 
\ 2 cos 2 i? 0 + 2 




(15) 


At. least up to i?o ~ 20° , the dependence T (i? 0 ) is insignificantly different from t he linear 
dependence T (i?o) = do/v , which simplifies the calculations. 

Since on the right side of Eq. (14) there are no terms that are linear in sin# , the interaction 
amplitude a (or the coil field value) can be any value other than zero, which affects only the time 
of the nutation damping rather than its character. 
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POSSIBLE REALIZATIONS 


To conclude the study of the nutation damping, let us make some final remarks relative to 
the application of the method proposed: 

a) As it is evident, the working part of the control moment of the n th coil pair 

~ sin (t/< + <r - /?„)■ At the same time, the projection of the rotor symmetry axis onto the 
direction r„ = z, cos ft n + z 2 sin /?„ , that corresponds to this pair is 

(s • r„) = T?cos(V' + O - 0n) + p COS ((7 - /?„) for small d and p . Since the moment transducer in 
the gyro incorporates four coil pairs, it is apparent that the control moment of any pair can he 
formed based on the signal taken from an adjacent pair. 

However, this mode of operation, in which the signal is taken from one pair and there are 
switchings-over in the other pair, can lead to an electrical engagement of both channels 
disrupting thereby their normal operation. In this respect, we need to obtain a signal with 
frequency t />* by anot her method, for example, using a. transducer to obtain a signal on one side 
of the rotor and a transducer to ensure control on the other side of the rotor. In this case, 
however, along with the moment effect, there will also be a force effect on the rotor or tension. 

b) Another method is the switching-over mode of operation, in which, for a certain period 
of time, the transducers are used for obtaining the information and then for control with 
periodic repetitions of the process. There is no input signal in the control mode, and it results in 
a decrease of K. Consequently, the control signal can be formed by a device of an underexcited 
oscillator type to the input of which an external signal with frequency V’* is initially applied. In 
this case, control will be efficient only over the time period until the phase difference bet ween 
the real signal from the rotor with frequency t/’* and the signal from the oscillator achieves the 
critical value (A<y? = tt/2) . 

If at the initial time it was the frequency , and the rotor deceleration could be simulated 
by the viscous friction when w(<) = u> 0 exp(-Al!) , the phase difference between the signals from 
the rotor and from the oscillator can be calculated as 

= J* u > o (l - e _Ar ) dr = w (t + je~ Xt - « ^ 0 A< 2 . 


If To = 27r/u> 0 , T x = A' 1 , then 

Tcon = ^7oT a (16) 

Let us consider the following example: if the rotor speed is U = 1000 r.p.s. 

(V>* = 2x • 1200 sec -1 ) , and the damping is such that 7\ = 1 hr , the nutation can be damped 
during T con « 1.22 sec until the phase difference becomes A tp = x/2 , after which we again need 
to switch over to the regime of signal reception and control moment formation. 

c) The proposed mechanism of nutation damping is active for both constant and variable 
magnetic fields generated by the coils, but the constant field efficiency is greater. The maximum 
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gain in efficiency of the effect on the rotor can be obtained, provided the transducer coils are 
used for control alone and the input signal with frequency ?/»* for the control moment formation 
is taken from other devices. 
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